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Fast algorithm computes histograms 

Lukasz Sliwczynski, University of Mining and Metallurgy, Cracow, Poland 



A histogram is a convenient 
method of presenting statistical 
features of data; you can use a his- 
togram for an estimate of the proba- 
bility density function of a random 
phenomenon. To prepare a histogram 
entails dividing the interval occupied 
by data values into some number of 
smaller intervals and then counting the 
number of occurrences of data in each 
subinterval. You can easily generate 
histograms with the aid of standard 
mathematical packages, such as Mat- 
lab. The problem lies in the fact that 
calculating a histogram using the stan- 
dard Matlab function hist{} from a 
large file of data (for example, 10 mil- 
lion data points) with high accuracy 
(for example, with 12-bit resolution) 
may take substantial time — many hours, 
or even days, depending on computer 
speed and memory. The problem is that 
hist{) treats its input data as 
"analog" data with almost infi- 
nite precision, and the routine must use 
program loops, which are very slow in 
Maflab. 

However, you often perform measure- 



LISTING 1 -MATLAB PROGRAM 
FOR FAST HISTOGRAM CALCULATION 



function N=fasthlst(g,L,W) 
% 

% N=fasthlst(g,L,W) 
% 
% 
% 
% 



Function calculates histogram of data from D/A converter. 
Range of numbers Is from to 2"L Subinterval length Is W 
L Is the number of bits of the converter. 



o=b(0'; 
gs=sort(g); 
gsd=dlff(gs); 
Index sflnd(gsd); 
v=[gs( lndex),gs(lengtli(gs))]+ 1; 
N=zfiros(l,2*L)i 

N(V)=[lodex(l),dlff(lndex),length(g5)-lndex(len 9 th(index))I; 
ffW~=l 

forl=i;flx(2 A L/W) 

N(l)=sum(N(W»(l-l)+ UW*Wt 

end; 

end; 

N=N(l:flx<2*L/W)>; 



ments with the use of A/D conversion, 
with the result that the input data for the 



Figure 1 



INPUT DATA (g) 

-r 



histogram-calculating 
procedure are discrete 
and that the smallest 
difference between any 
two data values is one. 
This fact allows you to 
write a special function 
for calculating the his- 
togram. You can use 
the complete Matlab 
program for calculat- 
ing histograms from 
discrete data (Listing 
1). The program is in- 
tended for placement 
into a Matlab m-file, 
which creates a func- 
tion from the routine. 
This function is acces- 
sible from the Matlab command window, 
just as any other function. Figure 1 illus- 
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Ordering of the data points in diagrams a through f follows the steps in Listing 1. 
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trates the operation of the algorithm. The 
first line ( 1 ) of the program prepares in- 
put data and gives the data the form of 
the column vector g (Figure la). The sec- 
ond line (2) performs data sorting and 
arranges the data in ascending order (gs 
in Figure lb). The difference vector gsd 
in Figure lc has peaks only at those 
points where data changes; elsewhere, the 
vector has a value of zero. The distance 
between two consecutive peaks is equal 
to the number of data samples of the 
same values, and the height of a peak is 
equal to the distance between two con- 
secutive numbers. 

Next, the routine (4) calculates the in- 
dexes of data changes (Figure Id) and 
calculates the values of data V (Figure le) 
appearing in input vector g (5). A"l" in 
this code line is necessary because a vec- 



tor index cannot equal zero in Matlab. In 
the following step (7), the step vector V 
indexes the output vector from the pro- 
cedure N. This indexing operation is sim- 
ple and fast because both indexes and 
data are natural numbers and belong to 
the same interval <0;2 L -1>. At this 
point, the essential part of the procedure 
is finished. The numbers N(i) (Figure If) 
are equal to the total number of occur- 
rences of data with the value i— 1 in the 
input vector g. The last six lines of the 
program reduce the initial resolution of 
the A/D converter that you use to acquire 
data. To obtain faster execution of the 
program, you switch on this part of the 
routine only when necessary — for exam- 
ple, when W does not equal 1. 

The effect of using this algorithm 
rather than the standard one in the Mat- 



lab package is significant. Using the stan- 
dard function hist{}, the time to compute 
a histogram from 10,000 points is 10 
sec; for the same data file with 1 2-bit i es- 
olution, the algorithm fasthist{} requires 
only 0.55 sec. These figures are based on 
using a 100-MHz Pentium machine with 
33 Mbytes of RAM, running under Win- 
dows 95. Using fasthist{}, a histogram 
from a 10 million-point data file t,-;kes 
approximately a half-hour as opposed to 
more than a day using hist{}. These re- 
sults occur using Matlab version 4.2. Ver- 
sion 5.2 gives somewhat different results: 
68 sec for hist{} and 0.38 sec for fasthist{}, 
a 179-times reduction. (DI #2422) 
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Loop powers current transmitter 

By N Kannan, Mediatronix Ltd, Triandrum, India 



In the circuit in Figure 1, the loop 
supply is 20 to 30V dc, and the loop 
current is 10 to 20 mA. IC, operates 
as a constant-current (L) source. IC 2 , 
Rj, and R 2 operate as a shunt regulator, 
which provides 12VV' . ICj is a TL032 
op amp. The sensor circuit could be a 
resistance- temperature-detector circuit 
or any other sensor- or signal-condi- 
tioning circuit. The IQ/Q, combination 
operates as a current sink, with 
I C =V IN /R E . The total loop cur 



rent is I p + I c - Because of its com- 
mon-mode input-range limitation, IC 2 
requires a minimum 2V input for prop- 
er operation. Hence, the V [N signal 
range must be 2 to 10V. With V |N = 2V, 
you adjust potentiometer P 1 to set I Loop 
at 10 mA. At this operating point, I c =4 
mA, and the balance of the current 
(I p =6 mA) goes to the shunt regulator, 
the op amp, and the sensor-circuit load. 
The maximum load current is typically 
5 mA. As V increases from 2 to 10V, 
I 100p increases proportionally, solely 
from the increase in I p , because I c is 
constant: 




12V 



Figure 1 
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Inexpensive components provide an efficient loop-powered current transmitter. 



rv- 



2V 

3V 2.4V 
4V 2.8V 
6V 3.6V 
8V 4.4V 
10V 5.2V 
For every IV increase in V JN , the 
change in V R is 0.4V. The circuit has good 



accuracy and temperature stability and 
operates with loop supplies of 2C to 30V. 
The maximum available sensor-circuit 
current is approximately 5 mA. Currents 
greater than this value result in poor V (r 
regulation. (DI#2421). 
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Spice creates time-variant resistor 

Vittorio Rkchiuti, Italtel, UAquila, Italy 



FIGURE 1 SHOWS A SUBCIRCUIT that 
represents the behavioral 
model of a linearly time- . 
variant resistor in PSpice (Oread, 
www.orcad.com). PSpice libraries in- 
clude no models of such resistors, which 
fou can use in various applications— for 
example, to simulate the fallback charac- 
teristic of a linear regulator in transient 
analysis. The suggested model uses the E- 
device in the PSpice analog-behavioral- 
modeling parts library. Listing 1 de- 
scribes the subcircuit in Figure I for a 
time-variant resistor with the character- 
istic in Figure 2. The function in Figure 
2 is the product of an independent volt- 
age source, V roNTROL . The independent 
zero-voltage source V E senses the cur- 
rent I(V SCs , SE ) through the resistor. The 
voltage V R between nodes 1 and 
2 of the time-variant resistor is 
V R : =f(t)*I(V SENSE ), where f(t) is a func- 
tion ofV CONTROL . (DI#2424) 



RESJJN 




EVALUE 
KV SENSE )*V(%IN+, %IN-) 



^CONTROL 



— 



A subcircuit provides the behavioral model of a time-variant resistor. 



is Figure 2 



CONTROL' 
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LISTING 1 -MACRO MODEL 
OF TIME- VARIANT RESISTOR 



SJBCKT 
R RI 



VAR„RES ! 2 
4 1G 



V_Vcomrol 4 
+FiVL2n Ik4n0 6n lk 

EJte-sJin 3 2 VALUE { I(V_Vsense)*V(4, 0) } 
* 

V_Vaim 1 3 DC OV AC QV ov 
.ENDS VAR RES 




The routine in Listing l simulates a linearly time-variant resistor. 



VCO supply touts low noise 

Ted Henderson, Linear Technology Corp, Milpitas, CA 

Many portable RF products use voltages that are higher than the primary herent disadvantages. The dc/dc convert- 

VCOs to generate the RF carrier battery voltage. Many designs use a dc/dc er produces noise that the regulator may 

frequency. These applications of- converter powering a low-noise linear not reject, resulting in regulator-output 

ten require low-noise VCO power-supply regulator. This approach has several in- noise far greater than the thermal-noise 
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re 1 



2.5 TO 
4.4V 



SHUTDOWN 
Q 



— j~ 4.7 p,F 2 



0.22 h,F;±: 



CPO 



^ 4.7 (iF 



VqUT 



SHDN 



LTC1682 fb 



GND 



4.2V 



36k- 



— 4.7 jiF 



r^lOOO pF 



»15k 




vco 

MURATA 
MQEOO1-902 



VCO 
OUTPUT 
902 MHz 



100k 
C — f-vW-O v c 



1000 pF 



1000 pF~ =p4.7nF 



I 



A single IC can generate a low-noise voltage for controlling a VCO, 



levels. The linear regulator might require 
a large output-compensation ca- 



Figure 2 



pacitor with stringent ESR re- 
quirements. Finally, the board area for 
both devices and support components 
can be large. The circuit in Figure 1 op- 
erates from an input range of 2.5 to 4.4V 
and generates a 4.2V low-noise voltage 
for a 900-MHz VCO. 

The circuit is based on IC,, an 
LTC1682 charge-pump dc/dc converter. 
The device includes a charge-pump/lin- 
ear-regulator tandem that's optimized 
for minimum regulator-output noise. 
The linear regulator operates with sev- 
eral types of output capacitors, includ- 
ing small, low-value, low-ESR ceramic 
capacitors. Figure 2 shows the close-in | 
phase noise of the VCO operating ^ 
in open -loop mode; Figure 3 
shows the peak-to-peak noise voltage at 
dbe regulator's output. (DI #2430). 
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AMPLITUDE 

(10 dB/DIV) 




CENTER = 902MHz 
SPAN = 100kHz 
SWP = 10s 
RESBW = 1kHz 
VBW = 30Hz 
REF = 0dBm 




The VCO in Figure l exhibits low close-in phase noise. 



(200 nV/DIV) 




100 i^SEC/DIV 

The supply in Figure l generates only approximately 200 mV of peak-to-peak noise. 
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IxC forms FM oscillator 

Abel Raynus, Armatron International, Melrose, MA 



A project required an inexpensive 
oscillator whose frequency increased 
step by step from 200 to 400 Hz and 
then decreased to 200 Hz. The first step 
was to design a VCO with a staircase driv- 
er. However, this approach entailed at 
least four ICs and many discrete compo- 
nents. An alternative method (Figure 1) 
requires only one 16-pin p,C (an 
MC68HC705KJ1, costing less than $1) 
and only a few external components. The 
process takes place exclusively in software 
(Listing 1). The u.C generates a burst of 
256 cycles of a given period. After that, 
the period decrements or increments, 
and the |xC generates a burst with the 
new period. The choice of decrement or 
increment depends on the contents of the 
flag register. If the flag is zero, it is a sign 
to decrement; otherwise, to increment. 
When the set of frequencies reaches com- 
pletion, the content of the register inverts. 
So, in one operation, the frequencies 
grow; in the next, they decrease. 

Keep in mind that a linear change in 
the value of the period causes a nonlin- 
ear change in the frequency value (ac- 
cording to a hyperbolic curve). So, you 
should perform all calculations in the 




Ci 
1.5 u.F 



I 

f"C3 4 MHz 
I I 



RESET A7 

MC68HC705KJ1 
0SC1 

0SC2 



ftl 



a OUTPUT 



An inexpensive replaces complicated VCO 
cirwrtry to configure a simple frequency-step 
generator. 



time domain, rather than in the frequen- 
cy domain. In our application, the non- 
linear aspect of the frequency is unim- 
portant. However, if necessary, you can 
realize any kind of FM by organizing the 
table of frequencies corresponding to 
each step. The design determines the pe- 



riod resolution, AT. In this case, it equals 
0.02 msec, because this value allows you 
to use only one 8-bit register (HPER) to 
perform all period-changing operations 
For more precision, you could use two 
or more registers. The number of fre- 
quencies in the set is K2 MAX =(T MAX — 



LISTING 1-fxC ASSEMBLY CODE FOR FREQUENCY GENERATOR 






1 


* Quasy FM Oscillator 








2 


* range 200 Hz - 400 Hz, 








3 


* period discretion - 0.02 ms 






4 


* number of frequencies - 


- 125 






5 


************* 


************************************* 






6 


*nolist 






0000 




7 


$include "std 


-j la. asm" 








6 


*list 










9 


•I/O PORTS 






0000 




10 


out equ 


7;prtA pinll 






11 


•CONSTANTS 






0000 




12 


HPmax equ 


250T; frequency i 200 Hz 






13 


•VARIABLES 






ooco 




14 


org 


RAM 




ooco 




15 


kl rmb 


1 ; sound duration counter 


00C1 




16 


MS rmb 


1 ; frequency number counter 


0OC2 




17 


HPER rmb 


1 ;Half 


of Period register 


00C3 




18 


flag rmb 


1 ,-flag register 






19 


•INITIALIZATION 




07F1 




20 


org 


MOR 




07F1 


20 


21 


feb 


900100000 ;Osc.parall.resitor 


0300 




22 


org 


ROM 




0300 


A6FF 


23 


init Ida 


t*f£ 




0302 


B704 


24 


sta 


ddrA ; set prtA for output 


0304 


3FC0 


25 


clr 


kl ;0 -> 


kl 


0306 


3FC1 


26 


clr 


k2 ;0 -> 


k2 


0306 


3FC3 


27 


clr 


flag 




03 OA A6FA 


28 


Ida 


HHPmax 




030C 


B7C2 


29 


sta 


HPER ;set T = 5 ms, f * 200 kHz 


O30E 


CD032D 


30 


main jsr 


burst; generate 256 cycles 


0311 


3DC3 


31 


tst 


flag 




0313 


2714 


32 


beq 


ml 




0315 


3CC2 


33 


inc 


HPER ; HPER + 1 


0317 


3CC1 


34 


m2 inc 


k2 ; k2 + 1 


0319 


B6C1 


35 


Ida 


k2 ; k2 - 


■> ACC 


031B 


A17D 


36 


crop 


8125T,- k2 < 125? 


031D 


25EF 


37 


bio 


main 




031F 


B6C3 


38 


Ida 


flag ; inverting flag 


0321 


A801 


39 


eor 


#i ; 




0323 


B7C3 


40 


sta 


flag ; 




0325 


3FC1 


41 


clr 


k2 ; -> k2 


0327 


20E5 


42 


bra 


main 




0329 


3AC2 


43 


ml dec 


HPER ; HPER 


- 1 


032B 


20EA 


44 


bra 


m2 




032D 


1E00 


45 


burst bset 


out,prtA 


1 -> out 


032P 


BEC2 


46 


ldx 


HPER 


HPER -> x 


0331 


CD0342 


47 


jsr 


dlyOlx 


half period delay 


0334 


1F00 


48 


bclr 


out, prtA 


-> out 


0336 


BEC2 


49 


ldx 


HPER 


HPER -> x 


0338 


CD0342 


50 


jsr 


dlyOlx 


half period delay 


033B 


3CC0 


51 


inc 


kl 


kl + 1 


033D 


3DC0 


52 


tst 


kl 


kl » 0? 


033F 


26EC 


53 


bne 


burst 




0341 


81 


54 


rts 




return from burst 


0342 


A602 


55 


dlyOlx Ida 


#2 ; delay O.Olx ms 


0344 


4A 


56 


repO decA 






0345 


26FD 


57 


bne 


repO 




0347 


5A 


58 


decx 






0346 


26F8 


59 


bne 


dlyOlx 




03 4A 


81 


60 


rts 




return from dlyOlx 






61 


************* 


************************************ 


07 FE 




62 


org 


VECTORS+6 




07FE 


0300 


63 


fdb 


init 
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T M)N )/AT. You can determine any fre- 
quency is the set as: 

FIag=0;f K =l/(T MAX -AT-K2). 

Flag=l; f K =l/(T MIN + AT-K2). 

Table 1 shows some of the frequencies. 
The duration of the set is t = 
^(T MAX +T, nN )-Kl MAX -K2 MfW . In this 
case, t = 2 msec (DI #2432). 
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TABLE 


20 
























nag 


K2 


1 (msec) 




HPER 


F(Hz) 


Comments 





1 


4.98 


1.49 


249 


200.8 







2 


4.96 


2.48 


248 


201.6 







3 


4.94 


2.47 


247 


202,4 











124 


2.52 


1.26 


12S 


396.8 










2.50 


1.25 


125 


400.0 





2.52 



1 3 



1.26 



126 



J9&.8 
393.7 
390.6 



Circuit converts pulse width to voltage 

Ron Hegler, Picker International, Cleveland, OH 



T 



he circuit in Figure I stems 
from a radio-con- 



trolled modeling ap- 



Figure 1 



plication, which requires a voltage 
proportional to the width of the 
incoming servo pulses. The circuit 
is optimized for a positive-going 
pulse width of 1 to 2 msec, re- 
peating at intervals of approxi- 
mately 17 msec. The output pro- 
duces a voltage of 0.95V for a 
1-msec pulse to 2.25V for a 2- 
msec pulse. The circuit operates 
similarly to a PLL, but it locks 
onto the pulse width, rather than 
to the frequency, of the incoming 
signal. IC JR is a one-shot multivi- 
brator with its time constant a 
function of R , the FET's on-re- 
sistance, and C., IC is a pulse- 
width comparator that compares 
the reference pulse from IC with 
the incoming pulse. Upon the rise 
of the incoming pulse, IC 's Q 
output clocks high and drives 
IC |B 's D input high. If LC times 
out before the input pulse falls, IC 1A 's Q 
output goes low, driving the D input of 
iC low. This action drives the Q of IC B 
high when the input pulse falls. The Q 




74HCQ4 



Vcc 





1- TO 2-mSEC TTL PULSE 
REPEATING AT 
APPROXIMATELY 60 Hz 



PR 

D Q 
ICia 
74HC74 



Ri 
51k 



Vcc 
O 



* AAA-, 



&10. 



PR 



D ICib Q 
74HC74 



CL Q 

Ol3 



— o 



Qi 

s 




R 2 

100k 



MPF102 



o.i Vf 



■ 1N4148 



O VOLTAGE 
OUTPUT 
1 mSEC=0,95V. 
-LT C 2 2 mSEC=2.25V. 
10 |xF 



This pulse-width-locked loop provides an output voltage proportional to the incoming pulse width. 



output connects to the FET through D, 
and R r C 2 filters the Q output to adjust 
the time constant of one-shot IC, , to 
match the incoming pulse. The voltage 



across C 2 indicates the incoming pulse 
width. (DI #2431). 
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Restore dc to NRZ sigmals 

Jay Kirschenbaum, Columbia University, New York, NY 

F 



or NRZ signals, ac amplification is I 
preferable to dc amplifica- 
tion, because ac amplifica- 



Figure I 



tion is usually more economical and has 
greater immunity to drift. However, ac 
processing has the disadvantage that ac 
coupling removes the dc reference level 
of the digital data stream. As Figure 1 
shows, if a long stream of ones or zeros 
appears at the amplifier's input, the out- 
put droops after a period determned by 
the RC time constant at the input. This 
problem makes it necessary to include a 
dc-restoration circuit following the ac 
amplifier. 

The circuit in Figure 2 restores the dc 
level that the amplifier removes and elim- 
inates droop, even if long strings of ones 



LONG STRING 
OF ONES 

n i — LrLP~U — — Li — ~LTiS~ 



> 



OUTPUT 
SIGNAL DROOPS 




Rom 



A long string of ones or zeros leads to droop in a purely ac -coupled digital amplifier. 



or zeros appear in the data stream. The 
circuit also provides additional amplifi- 
cation, increasing the signal to ECL lev- 



els. Synchronous clocking gives the cir- 
cuit good immunity to noise in the data 
channel. Tests reveal that the circuit op- 



Figurc 2 



JT-jm-JTJTTL 



(DATA STREAM 
FROM AC 
AMPLIFIER) 



.JTJTJTJTJTJTTL- 

CLK IN o- 

(CLOCK 
SYNCHRONOUS a. 
WITH DATA 50 
STREAM) 




C, 

0.047 nF 



JUlTLnilJl 



-OQout 



X 0.1 n-F 



NOTES: 

V EE =-5.2V. 
V +t =-2.0V. 

V D IS ADJUSTED AROUND - 1 ,3V. 
fig DEPENDS ON THE AMPLITUDE 
OF D, N . 



This circuit eliminates the droop from the data stream in Figure l and simultaneously amplifies the digital signal to ECL levels. 
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erates satisfactorily at clock speeds from 
dc to 1 GHz. 

IC is an MC10E416 five-channel line 
receiver, which boosts digital signals from 
a minimum of 50 mV p-p to ECL levels. 
The input circuit of each channel of ICj 
is a differential amplifier, which responds 
when the signal level at the true input 
swings past the signal at the comple- 
mentary input. In Figure 2, the input 
data stream goes to the_D channel. The 
complementary input D Q connects to the 
V M pin, which is at — 1.3V. The true in- 
put D connects to bias voltage V D 
through 5011. V D is also nominally 
- 1,3 V, but you can tweak it to compen- 
sate for device characteristics. IC 2 is an 
MC10EL31 flip-flop. Both ICs are speci- 
fied to 2 GHz. The logic level at D 1N trans- 
fers to Q OU[ at the positive transition of 
the signal at CLK IN . Resistor R } and the 
50fl pull-down resistor, R , form a volt- 
age divider that puts D above or below 
A/ D by half the expected input-voltage 
swing. The value of R } is a function of the 
swing of the input signal. For example , 



assume D has the value 100 mV p-p. If 
R } is 370X1, and Q OUT is at logic zero 
(- 1.75V), tiie voltage divider sets D Q 50 
mV below D , and Q 0UT stays at logic zero 
as long as these conditions prevail, 

Now suppose that logic one appears at 
Dw- This condition simultaneously puts 
D Q 50 mV above D , and Q oul switches to 
logic one { — 0.9V) at the next clock tran- 
sition. When this condition happens, the 
Rj-Rj voltage divider holds D Q 50 mV 
above D , and Q OUT stays at logic one un- 
til the next transition to logic zero. Volt- 
age droop disappears, and the dc refer- 
ence of the digital signal re-establishes 
itself. Moreover, the circuit amplifies the 
100-mV swing of the digital signal to 
ECL levels. Because the circuit triggers 
only when the clock edge appears, it's 
highly immune to noise in the data chan- 
nel. At power-up, one transition from 
logic zero to logic one and one transition 
from logic one to logic zero are necessary 
at the amplfier input to put the circuit in 
a defined state. After that sequence, the 
circuit faithfully follows the input. 



Proper timing between the data and 
the clock is important to reach the max- 
imum operating speed of the circuit. Ide- 
ally, the clock edge should arrive at ICj 
as soon as the data-setup time is over, 
switching Q OUT and pulling up Q as 
close to the data transition as possible. 
The theoretical upper limit on the oper- 
ating (clock) frequency is a function of 
the propagation delay through the circuit 
(825 psec at 25°C). Data transitions must 
be spaced at least this far apart, so the 
maximum clock frequency is 1.2 GHz. 
Tests show that the circuit operates suc- 
cessfully at clock frequencies of dc to 1 
GHz and data frequencies of dc to 500 
MHz. The tests used a square wave for the 
clock input and a synchronous square 
wave of a lower frequency at the data in- 
put. The circuit restores the dc level and 
amplifies the signal to ECL levels over the 
entire test range. (DI #2434). 

To Vote For This Design, 
Circle No. 364 
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TC 33370 Serie s A ppli 

• Consumer Electronics 

• Office Automation 

• Personal Computers 

« Wireless Communications 
- Industrial Systems 
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TC 33370 Series Featur es 

• Simplified Interface allows On/Off Control by a 
Microcontroller or other Circuitry 

• Lower System Cost with Reduced Component 
Count for a Total Converter System Solution 

• Improved Reliability with Internal Thermal 
Shutdown Protection, Cycle-by-Cycle Current 
Limiting, and Input Undervoltage Lockout with 
Hysteresis 

• Programmable State Controller 
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